Effects of anharmonic vibration on large-angle quasi-elastic scattering of ^^0+^''^Sm 
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We study the effects of double octupole and quadrupole phonon excitations in the ^'''*Sm nucleus 
on quasi-elastic '^®0-|-^'*''Sm scattering at backward angles. To this end, we use the coupled-channels 
framework, taking into account explicitly the anharmonicities of the vibrations. We use the same 
coupling scheme as that previously employed to explain the experimental data of sub-barrier fusion 
cross sections for the same system. We show that the experimental data for the quasi-elastic 
cross sections are well reproduced in this way, although the quasi-elastic barrier distribution has 
a distinct high energy peak which is somewhat smeared in the experimental barrier distribution. 
We also discuss the effects of proton transfer on the quasi-elastic barrier distribution. Our study 
indicates that the fusion and quasi-elastic barrier distributions for this system cannot be accounted 
for simultaneously with the standard coupled-channels approach. 

PACS numbers: 24.10Eq, 25.60.Pj, 25.70Bc, 27.60.-|-j 



I. INTRODUCTION 



The effect of channel coupling, that is, couplings of the 
relative motion between the colliding nuclei to their in- 
trinsic motions as well as transfer reactions, have been 
well known in heavy-ion collisions around the Coulomb 
barrier. In heavy-ion fusion reactions at sub-barrier en- 
ergies, the channel coupling effects enhance consider- 
ably the fusion cross sections as compared to the pre- 
diction of potential model calculation P, [3, It has 
been well established by now that the channel coupling 
gives rise to a distribution of potential barriers 0, [a, Q . 
Based on this idea, a method was proposed to extract 
barrier distributions directly from experimental fusion 
excitation functions by taking the second derivative of 
the product of center mass energy, and the fusion 
cross section, a[us{E), with respect to Coupled- 
channels calculations as well as high precision fusion 
data have shown that the fusion barrier distributions, 
^fus ^ (p[Ea{asiE)]/dE'^, is sensitive to the details of 
channel couplings, while the sensitivity is much more dif- 
ficult to see in the fusion cross sections [1, H, . 

Similar information as the fusion cross section can also 
be obtained from the quasi-elastic scattering (a sum of 
elastic, inelastic and transfer processes) at backward an- 
gles [10]. Timmers et al. measured the quasi-elastic scat- 
tering cross section for several systems 11[ , for which the 
fusion barrier distribution had already been extracted 
[^. They proposed that the corresponding barrier dis- 
tribution can be extracted by taking the first derivative 
of the ratio of the quasi-elastic to the Rutherford cross 
sections, daqci/daj^, with respect to the energy, E, i.e., 
f^"' = —d{daqci/daji)/dE. The properties of the quasi- 
elastic barrier distributions have been studied in more 
details in Ref [l^. These studies show that the quasi- 
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FIG. 1: Comparison between the experimental fusion (the 
filled circles) and quasi-elastic (the open squares) barrier dis- 
tributions for the ^^0-|-^'''*Sm reaction. They are normalized 
to unit area in the energy interval between -Bc.m. = 56 and 
70 MeV. The experimental data are taken from Refs. 0] and 

El. 



elastic barrier distribution is similar to the fusion barrier 
distribution, although the former is somewhat smeared 
and less sensitive to the nuclear structure effects. 

One of the systems which Timmers et al. measured is 
^^0-1- ^^''Sm [ll|. Figure 1 shows the comparison of the 
experimental barrier distribution extracted from fusion 
(the filled circles) and quasi-elastic (the open squares) 
processes. In order to compare the two barrier distribu- 
tions, we scale them so that the energy integral between 
Ec.ni.— 56 and 70 MeV is unity. For energies below 62 
MeV, the two barrier distributions resemble each other. 
However, at higher energies, they behave rather differ- 
ently, although the overall width of the distributions is 
similar to each other. That is, the quasi-elastic barrier 
distribution decreases monotonically as a function of en- 
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ergy while the fusion barrier distribution exhibits a dis- 
tinct peak at energy around i?c.m. = 65 MeV. So far, no 
theoretical calculations have succeeded in explaining this 
difference. The coupled-channels calculations of Tim- 
mcrs et al. [ll[ with the computer code ECIS 13'], which 
took into account the one quadrupole, 2^^, and the one 
octupole, 3^, phonon excitations of ^'*^Sm, were unable 
to reproduce both the experimental data of the quasi- 
elastic cross sections and the quasi-elastic barrier distri- 
bution. The ECIS results for the ratio of quasi-elastic 
scattering to the Rutherford cross sections fall off more 
steeply than the experimental data, while the obtained 
barrier distribution has a secondary peak similar to the 
fusion barrier distribution. They argued that this failure 
is largely due to the residual excitations not included in 
the ECIS calculations, which they postulated to be trans- 
fer channels. Esbensen and Buck have also performed the 
coupled-channels calculations for this system taking into 
account the second order couplings 14]. However, they 
did not analyze the quasi-elastic barrier distribution. 

These previous coupled-channels calculations took into 
account only the single phonon excitations in ^''^Sm. On 
the other hand, Hagino et al. psi . fl6t have shown that 
the double anharmonic quadrupole and octupole phonon 
excitations play an important role in reproducing the 
experimental fusion barrier distribution for this system. 
However, its effect on the quasi-elastic scattering has not 
yet been clarified so far. The aim of this paper is then 
to study weather the double anharmonic vibrational ex- 
citations of the "'^'^'^Sm nucleus can explain the difference 
in the shape of barrier distribution between fusion and 
quasi-elastic. The role of proton transfer reactions in this 
system is also discussed. 

The paper is organized as follows. In the next section, 
we briefly explain the coupled-channels formalism which 
takes into account the anharmonicities of the vibrational 
excitations. We present the results of our calculations in 
Sec. IH. We then summarize the paper in Sec. IV. 

II. COUPLED-CHANNELS FORMALISM FOR 
ANHARMONIC VIBRATION 

In this section, we briefly describe the coupled-channels 
formalism which includes the effects of anharmonic exci- 



tations of the vibrational states. We follow the proce- 
dure of Refs. [H, [11], which was successfully applied to 
describe the experimental fusion cross sections as well 
as the fusion barrier distributions of ^^0-|-^'''*'^^^Sm sys- 
tems. The total Hamiltonian of the system is assumed 
to be 



H = -— V2 + Hvib + Koup(r,e) (1) 



where r is the coordinate of the relative motion between 
the target and the projectile nuclei, ^ is the reduced mass 
and ^ represents the internal vibrational degrees of free- 
dom of the target nucleus. i?vib describes the vibrational 
spectra in the target nucleus. 

The coupling between the relative motion and the in- 
trinsic motion of the target nucleus is described by the 
coupling potential Vi^oup in Eq.(l), which consists of the 
Coulomb and nuclear parts. Using the no-Coriolis (iso- 
centrifugal) approximation 0, [l3| , they are given as 



Koup(r,C) = Vfc(r,C) + Viv(r,C), (2) 



Here O20 and O30 are the excitation operators for the 
quadrupole and octupole vibrations, respectively, and Rt 
is the target radius. The effect of anharmonicities for 
the quadrupole and octupole vibrations are taken into 
account based on the U(5) limit of the Interacting Boson 
Model (IBM). The matrix elements of the operator O — 
O20 + O30 in Eqs.© and O then read 
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for 6 low-lying states (i,j=l-6), where |1) = 10"*"), |2) — 
|2+), |3)_= |3-), |4) = |2+®2+), |5) = |2+ ® 3"), and 
|6) = |3 §53 ). In Eq.(5), (ii and /Ja are the quadrupole 
and the octupole deformation parameters, respectively, 
which can be estimated from the electric transition prob- 
abilities. The scaling of coupling strength with -x/iV, -^V 
being the number of boson in the system, is introduced to 
ensure the equivalence between the IBM and the geomet- 
ric model in the large N limit [18]. When all the x V^- 
rameters in Eq.(5) are set to be zero then the quadrupole 
moment of all the states vanishes, and one obtains the 
harmonics limit in the large IS! limit. Nonzero values of 
X generate the quadrupole moments, and, together with 
finite boson number, they are responsible for the anhar- 
monicities in the vibrational excitations. 



III. 



i6Q_^i44g]yj REACTION : COMPARISON 
WITH EXPERIMENTAL DATA 



We now apply the formalism to analyze the quasi- 
elastic scattering data of i6Q_|_i44g \{^ , The calcu- 
lations are performed with a version [l^ of the coupled- 
channels code CCFULL once the coupling matrix ele- 
ments are determined from Eq.(5). Notice that the iso- 
centrifugal approximation employed in this code works 
well for quasi-elastic scattering at backward angles p^ . 
In the code, the regular boundary condition is imposed 
at the origin instead of the incoming wave boundary con- 
dition. 



A. Effect of anharmonicities of nuclear vibrations 

In the calculations presented below, we include only 
the excitations in the ^"^^Sm nucleus whilst the excita- 
tions of the ^^O is not explicitly included. For sub-barrier 
fusion reactions, the latter has been shown to lead only 
to a shift of the fusion barrier distribution in energy with- 
out significantly altering its shape [13] , and hence can be 
incorporated in the choice of the bare potential. This is a 
general feature for reactions with the ^^O as a projectile. 
We have confirmed that it is the case also for the quasi- 
elastic barrier distribution. That is, although the ^^O ex- 
citations contribute to the absolute value of quasi-elastic 
cross sections themselves, the shape of quasi-elastic bar- 
rier distribution is not altered much. Since we are inter- 
ested mainly in the difference of the shape between the 
fusion and the quasi-elastic barrier distributions, we sim- 
ply do not include the ^^O excitations and instead adjust 
the inter-nuclear potential. 

For simplicity, we take the eigenvalues of the TJvib in 
Eq.(l) to be e = n^ti -\- 71363, where 712 and 713 are 
the number of quadrupole and octupole phonons, re- 
spectively. £2 and £3 are the excitation energies of the 
quadrupole and the octupole phonon states of the target 
nucleus, i.e., £2 — 1-61 MeV and £3 — 1.81 MeV, re- 
spectively. Notice that we assume the harmonic spectra 



I 10-^ 

D" 

D 
T3 



10"" 

0.20 

"0.15 
> 

(D 

."0.10 
0.05 
0.00 



CD 

cr 
Q 



1 1 1 1 

: (a) ^ 


1 1 1 1 

AHV(2ph) 

H0(2ph) ; 

HO(lph) 

^xp.(qel) 1 — • — 1 ■ 

^\ 

\ 'V 
\ 




"* 

\ \^ 
: \ m 

\ \ - 

\ \ 
\ \ 
\ \ 

1 \ 1 1 ' 


\ \ \ h 

- (b) 








7 \ W 




f \ T - 

/ ^V' 


\ 

\ 


/ V 


\ 

\ 

• V-... 







54 56 58 60 62 64 66 68 70 72 
Ec.m. [MeV] 



FIG. 2: Comparison of the experimental data with the 
coupled-channels calculations for ^'^O-f^^'^Sm reaction for (a) 
the ratio of quasi-elastic to the Rutherford cross sections and 
for (b) quasi-elastic barrier distribution. The dotted and 
dashed lines are obtained by including up to the single and 
the double phonon excitations in the harmonic limit, respec- 
tively. The solid line is the result of the coupled-channels 
calculations with the double anharmonic phonon excitations. 
The experimental data are taken from Ref. [lH . 



for the phonon excitations. It has been shown in Refs. 

[1^ that the effect of anharmonicity with respect to 
the excitation energy on the barrier distribution is in- 
significant once the energy of the single phonon states is 
fixed. The radius and diffuseness parameters of the real 
part of the nuclear potential are taken to be the same as 
those in Ref. 



i.e., ro = 1.1 fm and a = 0.75 fm, re- 
spectively, while the depth parameter is slightly adjusted 
in order to reproduce the experimental quasi-elastic cross 
sections. The optimum value is obtained as Vq = 112 
MeV. As usually done, we use a short-range imaginary 
potential with Wq = 30 MeV, r^, = 1.0 fm and — 0.3 
fm to simulate the compound nucleus formation. Finally, 

1 /3 

the target radius is taken to be Rt = 1.06 Aj( . We use 
the same values for the parameters /32, Ps, N,X2,X2fi and 
X3 as in Ref. [l5|. All the calculations presented below 
are performed at 9c.m. = 170°. 
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FIG. 3: (a) Comparison of the measured pure elastic (the 
open squares), the Z — 8 (—el) (the open circles) and the 
residual (the filled circles) components of daqci/dan with the 
coupled-channels calculations for ^''O+^^'^Sm reaction. The 
Z = 8 (—el) component is defined as the Z — 8 yields sub- 
tracted the elastic component, while the residual component 
the sum of Z = 6 and 7 yields. The dashed line is the result 
of elastic scattering, while the dotted line shows the inelastic 
cross sections for the single 2+ and 3~ phonon states. The 
solid line is the result of the sum of inelastic cross sections for 
the double phonon states in ^**Sm. (b) The same as (a) but 
for the pure elastic and the total inelastic cross sections. The 
experimental data are taken from Ref. pTj . 



The results of the coupled-channels calculations are 
compared with the experimental data in Fig. 2. Fig- 
ures 2(a) and 2(b) show the ratio of the quasi-elastic to 
the Rutherford cross sections, daqei/dan, and the quasi- 
elastic barrier distributions, D'^'^^ respectively. The dot- 
ted line denotes the result in the harmonic limit, where 
coupling to the quadrupole and octupole vibrations in 
^^^Sm are truncated at the single phonon level, i.e., only 
the 2+ and 3~ states are taken into account and all the 
X parameters in Eq.(5) are set to be zero. As we see this 
calculation fails to reproduce the experimental data. The 
obtained quasi-elastic cross sections, daqci/dcrji, drop 
much faster than the experimental data at high ener- 
gies. Also the quasi-elastic barrier distribution, _Di°', 



exhibits a distinct peak at energy around Ec.m. = 65 
MeV. These results are similar to the one achieved in 
Ref. [Tl|. The dashed line represents the result when 
the coupling to the quadrupole and octupole vibrations 
of ^"'^Sm is truncated at the double phonon states in the 
harmonic limit. In this case, we take into account the 
couplings to the 2+, 3", 2+ (g) 2+,2+ (g) 3" and 3" (g) 3" 
states. It is obvious that the results are inconsistent with 
the experimental data. To see the effect of anharmonic- 
ities of the vibrations, we then perform the the same 
calculations using the coupling matrix elements given in 
Eq.(5). The resultant quasi-elastic excitation function 
and the quasi-elastic barrier distribution are shown by 
the solid line. The calculated ratio of quasi-elastic to 
Rutherford cross sections quite well agree with the ex- 
perimental data. This suggests that the inclusion of an- 
harmonic effects in the vibrational motions is important 
for the description of the quasi-elastic excitation func- 
tions for the ^^0-|-^'*^Sm reaction. On the other hand, 
the result for D'^°^ is still similar to the barrier distribu- 
tion obtained by assuming the harmonic limit truncated 
at the one phonon level (the dotted line), although the 
former has a more smooth peak. 

Figure 3 shows the decomposition of the quasi-elastic 
cross sections to each channel for the calculation with 
the coupling to the double anharmonic vibrations (the 
solid line in Fig. 2). The fraction of cross section 
for each channel i in the quasi-elastic cross section, 
d<Ji/d<jqc\ — dui/\^^ rfo'i]; is also shown in Fig. 4. The 
open squares are the experimental elastic cross section 
while the open circles are the measured excitation func- 
tion for Z = 8 subtracted the contribution from the elas- 
tic channel. The latter contains not only the neutron 
transfer components but also the contributions of inelas- 
tic cross sections. The filled circles are the experimental 
residual (a sum of Z = 7 and Z = 6 yields) compo- 
nents of the daqc\/daR. The dashed line shows results of 
the coupled-channels calculations for the elastic channel. 
It reproduces reasonably well the experimental data for 
elastic scattering. The Z = 8 component of quasi-elastic 
cross sections is almost exhausted by the single phonon 
excitations, that is, the combined 2"*" and 3~ channels, 
as shown by the dotted-line in Figs. 3(a) and 4(a). The 
cross sections for the double phonon channels are given 
by the solid line in Figs. 3(a) and 4(a). These are im- 
portant at energies higher than around 66 MeV. If the 
components of all the inelastic channels included in the 
calculation are summed up, we obtain the dot-dashed line 
in Figs. 3 (b) and 4(b). 



B. Effects of proton transfer reactions 

In the previous subsection we have shown that the 
experimental quasi-elastic cross sections can be well ex- 
plained within the present coupled-channels calculations, 
which takes into account only the inelastic excitations in 
^'*'*Sm. However, the shape of quasi-elastic barrier distri- 
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FIG. 4: Same as Fig. 3, but for the fraction of cross section 
for each channel in the quasi-elastic cross sections. 

bution is still somewhat inconsistent with the experimen- 
tal data. As one sees in Figs. 3(a) and 4(b), the exper- 
imental data indicate that the charged particle transfer 
reactions may also play some role (see the filled circles in 
the figures). In this subsection, we therefore investigate 
the effects of proton transfer reactions, in addition to the 
anharmonic double phonon excitations. To this end, we 
use the macroscopic form factor for the transfer coupling 

dV(r) 
ar 

where Ft^ is the coupling strength and V{r) is the real 
part of the nuclear potential. In this paper, we con- 
sider a single proton transfer as well as the direct proton 
pair transfer reactions, although the experimental Z—6 
component may also include the alpha-particle transfer 
channel. The corresponding optimum Q— values for the 
transfer between the ground states are Qopt(lp) — —1.79 
MeV and Qopt(2p) = 0.13 MeV, respectively. The cou- 
pling strength Ft,, in Eq.(6) is determined so that the ex- 
perimental transfer cross sections for each Z=6 and Z=7 
components [22l | are reproduced. The optimum values 
for Ftr are found to be 0.12 and 0.16 fm for the one and 



FIG. 5: Effect of proton transfers on the quasi-elastic scatter- 
ing cross sections (the upper panel) and on the quasi-elastic 
barrier distribution (the lower panel) for ^®0-|-^'*''Sm reac- 
tion. The solid line is the result of the coupled-channels cal- 
culations including the effect of double anharmonic vibrations 
only. The dashed line is obtained by including, in addition, 
the couplings to the proton transfer channels. The experi- 
mental data are taken from Ref. [ill ]. 

the two proton transfer channels, respectively. 

The effects of proton transfer reactions on the quasi- 
elastic scattering is illustrated in Fig. 5. The solid line 
represents the results of the calculations including only 
the coupling to the double anharmonic vibrations. The 
dashed line is obtained by taking the coupling to the 
proton transfer channels into account, in addition to the 
anharmonic vibration channels. The upper panel shows 
the quasi-elastic cross sections, while the lower panel the 
quasi-elastic barrier distribution. We observe from Fig. 
5(a) that the inclusion of proton transfer reactions overes- 
timates the experimental daqci/dcrji at energies between 
62 and 68 MeV. Also the higher peak in the quasi-elastic 
barrier distribution becomes more distinct and thus wors- 
ens as compared to the calculation without the transfer 
channels. 

Figure 6 shows the contribution of each channel to the 
quasi-elastic cross sections. The fraction of each contri- 
bution is also shown in Fig. 7. The open squares are 
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FIG. 6: Contribution of quasi-elastic cross sections from sev- 
eral channels. The solid and dashed line are the results of the 
coupled-channels calculations for the proton transfer and the 
elastic cross sections, respectively. The dotted line denotes 
the sum of total inelastic and proton transfer cross sections. 
The corresponding experimental data are shown by the filled 
circles, the open squares, and the open triangles, respectively, 
which are taken from Ref. uM- 
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FIG. 7: Same as Fig. 6, but for the fraction in the quasi- 
elastic cross sections. 

the experimental elastic cross sections, while the filled 
circles and the open triangles are the experimental pro- 
ton transfer cross sections and the sum of total inelastic 
and transfer cross sections, respectively. The coupled- 
channels calculations for the elastic cross sections are 
shown by the dashed-line. Although it reproduces the ex- 
perimental data below around 62 MeV, it overestimates 
the data at higher energies. The sum of the contributions 
from the total inelastic and the proton transfer chan- 
nels is denoted by the dotted line, which reproduces the 
experimental data reasonably well, although the proton 



FIG. 8: Comparison of the theoretical fusion barrier distri- 
bution (dashed line) with the quasi-elastic barrier distribu- 
tion (solid-line) obtained with different coupling schemes for 
^®0-(-^^^Sm system. Both functions are normalized to unit 
area in energy interval between 54 and 70 MeV. (a) The re- 
sults of the coupling to one phonon state of quadrupole and 
octupole excitations of ^'*''Sm in the harmonic oscillator limit, 
(b) The same as (a) but for the coupling up to double phonon 
states, (c) The result when the coupling to anharmonic vibra- 
tion of double quadrupole and octupole excitations in ^*'*Sm 
is taken into account. 
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transfer cross sections themselves are underestimated at 
energies larger than 60 MeV (the solid line). The over- 
estimation of the quasi-elastic cross section indicated in 
Fig. 5(a) is therefore largely due to the contribution of 
elastic channel. 

From this study, we conclude that the inclusion of the 
proton transfer reactions in the coupled-channels calcula- 
tions does not explain the difference of the shape between 
the fusion and quasi-elastic barrier distributions for the 
i6o+i44gjj^ system. 

C. Discussions 

We have argued that the presence of high energy shoul- 
der, instead of high energy peak, in the quasi-elastic 
barrier distribution for the scattering between ^^O and 
^^^Sm nuclei cannot be accounted for within the present 
coupled-channels calculations, which take into account 
the anharmonic double phonon excitations in ^^"^Sm as 
well as the proton transfer channels. Figure 8 com- 
pares the calculated fusion barrier distribution D^'^^ and 
the corresponding quasi-elastic barrier distribution D'^'^^ 
for several coupling schemes as shown in Fig. 2 in the 
coupled-channels calculations. The solid line shows the 
quasi-elastic barrier distribution while the dashed line is 
for the fusion barrier distribution. They are normalized 
so that the energy integral between 54 and 70 MeV is 
unity. Figures 8(a) and 8(b) are obtained by including 
the one phonon and the two phonon excitations in ^^"^Sm 
in the harmonic limit, respectively. Figure 8(c) is the re- 
sult of the double anharmonic vibration coupling. From 
these figures, it is evident that the theoretical fusion and 
quasi-elastic barrier distributions are always similar to 
each other within the same coupling scheme, although 
the latter is slightly more smeared due to the low-energy 
tail [l^. This would be the case even with the excita- 
tions in ^^O as well as neutron transfer channels, which 
are not included in the present coupled-channels calcula- 
tions. Therefore, it seems unlikely that the experimen- 
tal fusion and quasi-elastic barrier distributions can be 
explained simultaneously within the standard coupled- 



channels approach. 



IV. CONCLUSION 

We have studied the effects of double anharmonic vi- 
brations of the ^''^Sm nucleus on the large angle quasi- 
elastic scattering for ^^0-|-^'''*Sm system. We have shown 
that the experimental data for the quasi-elastic scatter- 
ing cross sections for this reaction can be reasonably well 
explained. However, we found that the obtained quasi- 
elastic barrier distribution still shows the clear doubled- 
peaked structure, that is not seen in the experimen- 
tal data. This was not resolved even if we took the 
proton transfer channels into account. Our coupled- 
channels calculations indicate that, within the same cou- 
pling scheme, the quasi-elastic and fusion barrier distri- 
butions are always similar to each other. Although de- 
tailed analyses including neutron transfer channels in a 
consistent manner are still necessary, it is thus unlikely 
that the fusion and quasi-elastic barrier distributions can 
be explained simultaneously with the standard coupled- 
channels framework. This fact might be related to the 
large diffuseness problem in sub-barrier fusion, in which 
dynamical effects such as couplings to deep-inelastic scat- 
tering are one of the promising origins [Ij, lU, . It is 
still an open problem to perform the coupled-channels 
calculations with such dynamical effects and explain the 
difference of the shape between the fusion and the quasi- 
elastic barrier distributions for the -'^^0-|-^'*^Sm reaction. 
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